Reaction cross sections of carbon isotopes incident on a proton 
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We systematically study total reaction cross sections of carbon isotopes with N — 6-16 on a proton 
target for wide range of incident energies, putting an emphasis on the difference from the case of 
a carbon target. The analysis includes the reaction cross sections of 19 . 2 °. 22 Q a t 40 AMeV, the 
data of which have recently been measured at RIKEN. The Glauber theory is used to calculate the 
reaction cross sections. To describe the intrinsic structure of the carbon isotopes, we use a Slater 
determinant generated from a phenomenological mean-field potential, and construct the density 
distributions. To go beyond the simple mean-field model, we adopt two types of dynamical models: 
One is a core+n model for odd-neutron nuclei, and the other is a core+n+n model for 16 C and 22 C. 
We propose empirical formulas which are useful in predicting unknown cross sections. 
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Reactions of unstable neutron-rich nuclei with a proton 
target are of current interest [l|, 0| since such reactions 
are at present the major means to sensitively probe the 
matter densities of exotic nuclei, especially the region 
of nuclear surface. If one appropriately selects incident 
energies, protons could be more sensitive to neutron dis- 
tributions than proton distributions of nuclei. 

The structure of carbon isotopes has recently attracted 
much attention. Several works have been done already 
experimentally 0, H i, H and theoretically @, & ©, HO, 
Hi] . For example, the structure of 22 C has been studied 
by two (W.H. and Y.S.) of the present authors in a three- 
body model of 20 C+n+n. They showed that it has a 
Borromean character [ijj • 

The purpose of this paper is to report a systematic 
analysis of the total reaction cross sections of carbon iso- 
topes incident on a proton at energies from 40 AMeV to 
800 AMeV, and predict the cross sections of neutron-rich 
isotopes. We also estimate the cross sections contributed 
by protons or neutrons in the nuclei of carbon-isotopes. 
This study is motivated by an ongoing measurement of 
the reaction cross section of 22 C at RIKEN (l3| . 

Recently, we have performed systematic analyses of to- 
tal reaction cross sections of carbon isotopes on 12 C for 
wide energy range using the Glauber model [l4]. We 
found reasonable parameterizations of nucleon-nucleon 
scattering amplitudes, and obtained fairly good agree- 
ment with available data. We predicted the total reac- 
tion cross section of a neutron-rich isotope, 22 C on 12 C, 
and obtained a sizable effect of the extended surface. An- 
other purpose of this study is to discuss the advantage 
and disadvantage of a proton and a carbon target. 

In this paper, we adopt the same prescription as our 
previous work [lj] for describing the nuclear structure, 
and calculate the total reaction cross sections of proton- 
carbon isotopes similarly to the case of 12 C target. Wc 



calculate systematically total reaction cross sections for 
wide energy range using the Glauber model. Of course, 
we should note that the model may not be so good at 40 
AMeV. 

We treat the interactions of proton-proton and proton- 
neutron separately. The wave functions of carbon iso- 
topes are generated based on a simple mean-field model. 
To go beyond that, we adopt two types of dynamical 
models: One is a core+n model for an odd N nucleus, 
and the other is a core+n+n model for 16 C and 22 C. The 
reason for the latter model is explained in Ref. [l4|. We 
do not take into account the Coulomb potential, which 
would affect the magnitude of the cross sections for the 
low energy processes to some extent, but, for the present 
discussion, the effect is minor. 

This paper is organized as follows: The reaction mod- 
els for the calculations of reaction cross sections are pre- 
sented in Sec. II. We explain our input data in Sec. III. 
We present the cross section calculation in Sec. IV. The 
contributions of the protons and neutrons inside an iso- 
tope to the reaction cross section are presented in Sec. V. 
Summary is given in Sec. VI. In Appendix, we discuss 
the parameterization of the nucleon-nucleon scattering 
amplitude. 



II. THE GLAUBER MODEL FOR REACTION 
CROSS SECTION CALCULATIONS 

Here we summarize our basic formula for the following 
discussions. 

The total reaction cross section of proton-nucleus col- 
lisions is expressed as 
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where b is the impact parameter vector perpendicular 
to the beam (z) direction, and x(^) is the phase-shift 
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function denned below. We calculate this quantity using 
the Glauber theory. 

The Glauber theory provides us with an excellent 
framework to describe high energy reactions. In this 
framework, the optical phase-shift function (the elastic 
S'-matrix) for proton-nucleus scattering is given by flEI ] 



ix(b) _ 



r pp (b + Sl ) |Vo>, (2) 



where is the intrinsic (translation-invariant) A- 
nucleon wave function of the projectile's ground state 
(A is the mass number of the projectile), and s, is the 
projection onto the a;?/-plane of the nucleon coordinate 
relative to the center-of-mass of the projectile. Here 
is 1 for neutron and —1 for proton. 

When we apply this framework to low energy pro- 
cesses, such as the one less than 100 MeV, its usefulness 
should be carefully assessed. As a prescription, we care- 
fully choose the parameters of the nucleon-nucleon scat- 
tering amplitude so as to reproduce the reaction cross sec- 
tions of proton- 12 C scatterings in consistent with those 
of 12 C- 12 C scatterings [lj]. 

The profile function, T p n, for pp and pn scatterings, is 
usually parameterized in the form; 



r P Ar(b) 



i 



ia P N 



Att(3 pN 



a P N e 



-b 2 /(20 pN ) 



(3) 



where a P N is the ratio of the real to the imaginary part of 
the pp {pn) scattering amplitude in the forward direction, 
a P °N ^ s t ne PP (p n ) total cross sections, and P p m is the 
slope parameter of the pp (pn) elastic scattering differen- 
tial cross section. We parameterize the nucleon-nucleon 
scattering amplitude with a single Gaussian, because we 
find that double Gaussians give numerically almost the 
same reaction cross sections as the single Gaussian. We 
discuss this point in Appendix. 

There are several approximate expressions of the 
Glauber model on the market. We explain some of the 
expressions below. 

In the optical limit approximation (OLA), the phase- 
shift function of proton-nucleus scattering is given by 



3 iXOLA(b) _ 
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with 



iX P (b) = -j drp p (r)T pp (s + b), 

iXn{b) = - j ' drp n (r)Tpn(s + b), (5) 

where Xp (Xn) implies the phase shift due to the pro- 
tons (neutrons) inside the nucleus. The function p p (r) is 
the proton density distribution, and p n (f) is the neutron 
density. 



In the few-body (FB) calculation, the OLA is used 
for the integration involving the coordinates of the core 
nucleons, while the integration for the valence- nucleon 
coordinate is performed without any approximation (l6l . 
Ell EE El- In this treatment, Eq. (2]) is reduced to the 
following expression for the case of core+n configuration: 
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with 
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where tpo is the single-particle wave function of the va- 
lence nucleon and be is the impact parameter between 
the proton and the core. The phase-shift function, xcp, of 
the proton-core scattering is defined in exactly the same 
way as Eq. Q. The proton-neutron phase-shift func- 
tion, Xpn, is defined through the relation; exp(ix pn (b)) = 
1 — r p „(6). In this paper, we adopt both Eqs. (fj| and 
© to calculate the phase-shift function. 

In the discussion below, we adopt the kinematics of 
the projectile's rest frame. We specify processes by the 
energy of an incident proton. For example, the energy of 
40 ^4MeV of an incident nucleus of the mass number A 
in the proton-fixed frame corresponds to the energy of 40 
MeV of an incident proton in the projectile's rest frame. 



III. INPUT DATA 

In this section, we list the input quantities in the cal- 
culations of the reaction cross sections, and give some 
discussions. 

The inputs for Eq. |2]) are the projectile's intrinsic A- 
nucleon wave function and the parameters of the pp and 
pn profile functions. For Eq. ([1]), we need only proton 
and neutron intrinsic densities and the parameters of the 
profile functions. For Eq. ©, we need the single-particle 
wave function for the valence neutron as well. 
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are taken from 



[20|| , the experimental values of 



Cpp, Vpni a pp an d a pn are listed in energy range from 20 
MeV to 300 MeV. 
The parameters (3 PP 



and P p 



j pn are determined from the 
fact that the total elastic cross section, cr pN , is equal to 
the total cross section in this energy range, since only the 
elastic scattering is energetically possible until the pion 
production threshold is open. For the profile function, 
Eq. we have ^ 



cl _ 1 + a pjV / totV 



(8) 



Since a^ N = a p °^ , we can derive the following expression 
for P pN ; 
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The experimental data of Ref. 20] and the f3 P N values de- 
termined from Eq. © are listed in Table [H In Ref. [2l| . 
all the needed parameters are listed in the energy range 
from 100 MeV to 1000 MeV. They are also given in Ta- 
bled] For the energy higher than 300 MeV, (3 p n is deter- 
mined from Eq. ((8|) using both the data of a^ N and offi. 
The data on a° l N are taken from PDG tabulation [22| 
and the uncertainty of the data is fairly large. 

As we show in Appendix, the description of the pn 
elastic differential cross section with these parameters is 
reasonable, but not perfect especially in the forward di- 
rection. Fortunately, this does not affect the total reac- 
tion cross sections. We discuss it in some detail in the 
appendix. 

The densities that we use here are constructed from a 
core+n model for the odd isotopes, 13 > 15 > 17 > 19 C, where the 
cores are 12 > 14 . 16 > 18 C, respectively. For 16 ' 22 C, a core+2n 
model is assumed. The densities of the carbon isotopes 
are displayed in Fig.[TJ and the corresponding root-mean- 
square (rms) radii are summarized in Tabic HT1 The detail 
of these densities can be found in Ref. [Tl |. 

IV. PREDICTION OF THE REACTION CROSS 
SECTIONS 

Here we show our numerical results of the total reac- 
tion cross sections of proton-carbon isotopes reactions. 

Before we predict the reaction cross sections for the 
isotopes, we first show how well our densities and the pa- 



TABLE I: Parameters of the pn and pp profile functions as 
defined in Eq. @. E is the projectile's incident energy. 



E 


tot 
V pp 






tot 
O pn 


(%pn 


ftpn 


(MeV) 


(fin 2 ) 




(fm 2 ) 


(fm 2 ) 




(fm 2 ) 


40 


7.0 


1.328 


0.385 


21.8 


0.493 


0.539 


60 


4.7 


1.626 


0.341 


13.6 


0.719 


0.410 


80 


3.69 


1.783 


0.307 


9.89 


0.864 


0.344 


100 


3.16 


1.808 


0.268 


7.87 


0.933 


0.293 


120 


2.85 


1.754 


0.231 


6.63 


0.94 


0.248 


140 


2.65 


1.644 


0.195 


5.82 


0.902 


0.210 


160 


2.52 


1.509 


0.164 


5.26 


0.856 


0.181 


180 


2.43 


1.365 


0.138 


4.85 


0.77 


0.154 


200 


2.36 


1.221 


0.117 


4.54 


0.701 


0.135 


240 


2.28 


0.944 


0.086 


4.13 


0.541 


0.106 


300 


2.42 


0.626 


0.067 


3.7 


0.326 


0.081 


425 


2.7 


0.47 


0.078 


3.32 


0.25 


0.0702 


550 


3.44 


0.32 


0.11 


3.5 


-0.24 0.0859 


650 


4.13 


0.16 


0.148 


3.74 


-0.35 


0.112 


700 


4.43 


0.1 


0.16 


3.77 


-0.38 


0.12 


800 


4.59 


0.06 


0.185 


3.88 


-0.2 


0.12 


1000 


4.63 


-0.09 


0.193 


3.88 


-0.46 


0.151 



rameters of the profile functions fit the experimental data 
of the proton- 12 C total reaction cross sections. In Fig. [21 
we compare the numerical results with the experimental 
data over the energy range from 40 MeV to 800 MeV. 
As one can see from the figure, they reasonably agree 
with each other over all the energy range. At energies 
lower than 100 MeV, where the data fluctuate by 15 % 
at most, our results follow the largest data. In this en- 
ergy region, a systematic uncertainty of our approach is 
estimated to be about 15 %, which is consistent with the 
estimation by two of us (B.A.-I. and Y.S.) for the case of 
6 He+ 12 C reaction 40 AMeV [Hj]. They confirmed that 
the eikonal approximation gives about 15 % larger cross 
sections than those by the quantum-mechanical (exact) 
calculation. 

Now we show our predictions for all the carbon isotopes 
at all the energies using the parameters given in Table [U 
The numerical results of the total reaction cross sections 
are summarized in Table [HI] 

Let us estimate the contributions of the breakup ef- 
fect although it is expected to be small for a proton tar- 
get. Equations @ and $Q contain the breakup effect, 
while Eq. ^ does not. We compare them to estimate 
the breakup effect. As an illustrative example, we calcu- 
late the reaction cross section of a typical halo nucleus, 
19 C, incident on a proton using Eqs. ((4|) and ((6]). We as- 
sume the structure of 19 C as 18 C+n with the one-neutron 
separation energy of 0.581 MeV [24j . The numerical re- 
sults at 40 MeV and 800 MeV are 763 mb (758 mb) and 
372 mb (373 mb) respectively when Eq. © (Eq. |@J) is 
used. The difference is less than one %, which is consis- 
tent with the results of Ref. [l6| . The breakup effect can 
therefore be neglected. This validates our use of Eq. 

For convenience, we introduce the black-sphere radius, 
a, defined through [25[ 

(T R = na 2 . (10) 

Following the Carlson's prescription |26j, we fit the nu- 
merical results by parameterizing the radius, a, using a 



TABLE II: The rms radii in fm of matter, neutron and proton 
density distributions for the carbon isotopes. 



Isotopes 


<m 


r n 


r p 


12 C 


2.31 


2.30 


2.33 


13 c 


2.37 


2.40 


2.34 


14 c 


2.39 


2.46 


2.31 


15 C 


2.65 


2.84 


2.34 


16 c 


2.66 


2.83 


2.34 


17 c 


2.94 


3.20 


2.38 


18 c 


2.78 


2.96 


2.36 


19 c 


3.09 


3.37 


2.38 


20 c 


2.99 


3.23 


2.37 


22 C 


3.58 


3.92 


2.43 
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r (fm) r (fm) 



FIG. 1: Density distributions of the carbon isotopes. The dotted curve shows the proton density, the dashed curve the neutron 
density, and the solid curve the matter density. 
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FIG. 2: Comparison of the numerical results with the experi- 
mental data for the total reaction cross sections of proton- 12 C 
reaction as a function of energy. The experimental data are 
taken from Refs. [HIH- 

TABLE III: Total reaction cross sections of proton-carbon 
isotopes in units of mb. E is the projectile's incident energy. 



E 
(MeV) 


12 C 


13 c 


14 c 


15 C 


Isotopes 

16 C 17 C 


18 C 


19 C 


20 C 


22 C 


40 


432 


467 


489 


580 


605 


682 


662 


758 


761 


957 


100 


284 


308 


327 


372 


394 


436 


443 


491 


509 


604 


200 


218 


236 


252 


282 


300 


330 


340 


372 


390 


453 


300 


202 


218 


231 


257 


273 


299 


309 


337 


353 


407 


425 


200 


214 


227 


251 


265 


289 


300 


324 


339 


389 


550 


217 


231 


242 


267 


282 


307 


315 


341 


356 


408 


650 


233 


247 


259 


284 


299 


324 


332 


359 


374 


429 


800 


243 


257 


268 


294 


309 


335 


342 


373 


385 


442 



simple geometric picture with a correction term; 

a = C + r A 1 ^. (11) 

This includes a A 1 / 3 correction in addition to the simple 
geometrical A 2 / 3 term. In Ref. [26| . Carlson used R p 
instead of Co here. He fitted the reaction cross sections 
of stable nuclei incident on a proton target in the energy 
range from 40 MeV to 560 MeV. 

In Fig. [3l we compare our numerical results (open cir- 
cles) listed in Table [01] with the fit using Eq. (JTTJ) (solid 
lines) at 40, 100 and 550 MeV. The values of C and 
ro extracted from the fit are given in Table IIVI These 
values that we find are different from those obtained by 
Carlson, which are given in the parentheses in the table. 

The parameter Co implies the strength of A 1 / 3 correc- 
tion to A 2 / 3 -dependence of the reaction cross sections. 
The values in Table IIVI decrease with the energy, which 
is consistent with the geometrical picture of the cross 
section at high energy, because ctr oc A 2 ! 3 for proton- 
nucleus reaction. 
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FIG. 3: Comparison of the black square radii, a = y ctr/tt, 
for the numerical results of the carbon isotopes (open circles) 
with the fit using Eq. (JTTJ (solid line) as a function of the 
mass number. 

As one can see from Fig. [31 the curves with our pa- 
rameters nicely fits the numerical results for the stable 
isotopes as well as the neutron-rich unstable isotopes. 
The estimations with Carlson's parameters underesti- 
mates our numerical results for neutron-rich cases. This 
reflects an anomalous mass number dependence of the 
size of such exotic nuclei. Even with the new parameter- 
ization, the reaction cross section of 22 C is even larger 
than the fit, especially at 40 MeV. This would suggest 
an extended surface structure of 22 C. We believe that 
this simple fitting formula will serve as a reference for 
discussions of the total reaction cross sections. 

Moreover, we empirically deduce the rms nuclear mat- 
ter radii, using the black-sphere radius, a, of Eq. (fTOj) . If 
we assume a rectangular density distribution for nuclei, 
we obtain 

This radius clearly depends on the incident energy. At 
around 100 MeV, this tbs value happens to agree rea- 
sonably well with r m listed in Table |TTJ This suggests 
that we may empirically access to the rms nuclear mat- 
ter radii of carbon isotopes just by measuring ctr at 100 
MeV. This is consistent with the estimations in Ref. [25| . 
The authors of Ref. [H[ pointed out that, for T p > 800 
MeV, tbs almost completely agrees with the empirically 
deduced values of the rms matter radius for stable nu- 
clei having mass A > 50, while it systematically deviates 
from the deduced values for A < 50 [25}. Since carbon 
isotopes belong to light nuclei, we may choose the energy 
which gives a little bit larger ctr, to obtain trs close to 

' 171 

We propose another empirical formula. For all the car- 
bon isotopes, we find that the following relation is satis- 
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fled over all the energy range: 



wv c) 



R(C) 



g tot 



N tot 
w pn 



5 a tot I Q^tot ' 
^ w pp 1 uu pn 



(13) 



with R(C) = 0.96 ± 0.05. Here N > 7 and (^(c^n) 
is the proton-proton (proton-neutron) total cross section 
at a given energy. The value of R(C) is obtained by av- 
eraging the 153 numerical results (9 isotopes times 17 
energy points) of the reaction cross sections, and 0.05 is 
the standard deviation of these points. At high energy, 
i?(C) of Eq. (fT3|) is very close to unity. In Fig. |4j we 
plot R(C) for selected energies. Only at 40 MeV, some 
points come slightly below this relation, which would sug- 
gest the breakdown of the approximations, such as the 
fixed-scatterer approximation, contained in the Glauber 
model. 

At least for carbon isotopes, the expression (JT3J) indi- 
cates that if we know the reaction cross section of a stable 
isotope, we can predict the reaction cross section of other 
isotope within the error bar. Whether this holds for any 
nuclides or not is left for a future study. 

Experimental data are expected to appear at around 40 
MeV for 19. 2 0.22 C [7J] Here we pre dict them. In Fig. 
we compare our prediction for the reaction cross sections 
of carbon isotopes with the available experimental data. 
The preliminary data of proton- 22 C reaction cross section 
has been re por ted to be around 1000 mb with a large 
uncertainty [131 ], which is consistent with our prediction. 

For 22 C, we generate several densities that give dif- 
ferent two-neutron separation energies of 0.489, 0.361, 
0.232 and 0.122 MeV for the last two neutrons. All of 
them lie within the error bar of the experimental value, 
0.423± 1.140 MeV [IJ]. Using these densities, we calcu- 
late the reaction cross sections for proton- 22 C reaction 
at 40 MeV in order to examine the separation-energy de- 
pendence. The results are listed in Table IVl The change 
in radius from 3.6 to 4.1 fm gives change in the reaction 



TABLE IV: The parameters of Eq. (JTTJ) which give the lines 
plotted in Fig. The values in the parentheses are those of 
Carlson [H. 



E (MeV) 


Co (fm) 


r (fm) 


40 


-3.83 (1.00) 


3.27 (1.21) 


80 


-3.123 


2.73 


100 


-2.95 (-0.31) 2.58 (1.37) 


140 


-2.68 


2.38 


200 


-2.46 


2.21 


240 


-2.36 


2.14 


300 


-2.14 


2.03 


425 


-1.62 


1.85 


550 


-1.58 (-0.30) 1.84 (1.33) 


800 


-1.31 


1.782 



cross section of about 50 mb at 40 MeV. At 800 MeV, the 
change in the reaction cross section is about 10 mb. This 
gives an estimate of an uncertainty of our calculations. 



V. THE NEUTRON CONTRIBUTION TO THE 
REACTION CROSS SECTION 

Here we estimate that contribution to reaction cross 
sections which comes from the neutrons in the nucleus. 
Reactions with a proton are superior to those with a 12 C 
when we look into such separate contributions, because 
12 C is equally sensitive to protons and neutrons. For the 
purpose of discussion here, we use the following relation: 



i _ | giXn (*>)+%>(*>) | 2 = | e 'Xp(f>)| 2 ^1 . 

+ le^WlVl 



l_| e W>)| 2 ) (l-| e ^ 6 ) 



(14) 

|2 X 



Then we define the proton-nucleus reaction probability 
Pa (b) and its decomposition into neutron and proton 
contributions, P n (b) and P p (b), as 



P A (b) = P n (b) + P p (b), 



(15) 



where 



iX P (b) 



C 1 



Pn(b) 



P p {b) = |e lx "^ 
+ d[l 



e ix»(i>) 
JXn(b) I 2 



ix»(&) 

1 - 
'X P (b)| 2 
1 - I 



iXp(b) 



JX P (b)\ 2 



(16) 



where c + d = 1 , and c and d represent the neutron and 
proton contributions from the interference term, respec- 
tively. Equation (|Tj) is expressed as 



ctr = 2tt / bdbP A {b). 
Jo 



(17) 



TABLE V: The reaction cross sections of 22 C incident on a 
proton target at 40 AMeV for different two-neutron separa- 
tion energies, S2n- The r m value denotes the rms matter 
radius. 



S 2n (MeV) 


r m (fm) 


ctr (mb) 


0.489 


3.6 


957 


0.361 


3.7 


969 


0.232 


3.8 


985 


0.122 


4.1 


1005 
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1.2 
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0.8- 
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800 MeV 



~l 1 1 1 1 

14 16 18 20 22 
A 



1.4- 
1.2 

1 — , 
0.8- 

0.6 



550 MeV 



~l 1 1 1 1 

14 16 18 20 22 
A 



1.4- 
1.2 

1-0—1 
0.8- 
0.6 



425 MeV 



□ □□□□□□ 



"1 1 1 1 1 

14 16 18 20 22 
A 



1.4- 
1.2 

1-0-1 
0.8- 
0.6 



100 MeV 



"1 1 1 1 1 

14 16 18 20 22 
A 



1.4- 

1.2- 

1.0 

0.8- 

0.6- 



40 MeV 



□ □ □ □ 

U -g P a. 



1 1 1 1 1 

14 16 18 20 22 
A 



FIG. 4: The values of R(C) of Eq. (|13p as a function of the mass number, A. The dashed lines denote 0.96±0.05 values. 



The values of c and d must satisfy the condition, 
c + d = 1, but the choice of them is not unique. Here 
we discuss two cases; 1) c = d = 1/2, 2) c = N/A and 
d = Z/A, to see the dependence of the choice. The for- 
mer implies that both neutrons and protons contribute 
equally to the interference term, while the latter implies 
that the contribution of the neutrons and protons to the 
interference term is proportional to their numbers. 

In Fig. [5J we show our predictions of the total re- 
action cross section of proton- 22 C reaction as a func- 
tion of energy (solid curve). We also draw the neu- 
tron and proton contributions for c = d = 1/2 (dash- 
dot-dotted and dash-dotted) and c = N/A, d = Z/A 
(dashed and dotted), respectively. Due to the fact that 
22 C is very neutron-rich, we learn from this figure that 
the neutron contribution dominates the reaction cross 
sections. Also we find that the proton and neutron con- 
tributions depend modestly on the choice of the values 
of c and d. For example, for c = d = 1/2, the neutron 



contribution to the total reaction cross section of 22 C, 
2tt r o °° bdbP n {b)/2TT r o °° bdbP A (b), is about 0.87 and 0.73 
at 40 MeV and 800 MeV, respectively. For c = N/A, 
d = Z/A, the neutron contribution to the total reaction 
cross sections is about 0.93 and 0.80 at 40 MeV and 800 
MeV, respectively. 



According to Eq. (jl3|) . the neutron contribution to the 
total reaction cross section would be similar to the ratio 
N(T™/(Zct$* + No$). At 40 MeV and 800 MeV, the 
ratio of Na^/iZa^f + Naffl reads 0.89 and 0.69, re- 
spectively. These values are quite similar to the above 
ratios of our numerical results. 

Figure [7] shows the reaction probability times 27t6 of 
proton- 12 C reaction as a function of the impact param- 
eter, b. We plot 2nbP(b), because this quantity more 
directly reflects the contribution to ctr, than P(b) itself. 
The solid curve represents the total reaction probability, 
Pa (b), in Eq. (fH))) . The neutron contribution P n {b) is 
shown by the dashed curve and the proton contribution 



8 





1000 -i 




900 - 




800 - 




700 - 


E_ 


cc 


600 - 


B 




500 - 




400 - 




300 - 



□ □ 



10 12 14 16 18 20 22 24 

A 

FIG. 5: Reaction cross section for the carbon isotopes at 40 
MeV. The experimental data are taken from Ref. [26]. The 
larger one is natural carbon at 42 MeV, and the smaller one 
is 12 C at 40 MeV. The energy is converted to the case of a 
proton target. The preliminary data for 22 C is about 1000 
mb with a large uncertainty [la ]. 
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FIG. 6: The total reaction cross section of proton- 22 C as a 
function of energy (solid curve), and its decomposition to the 
neutron and proton contributions. 



P p (b) is shown by the dotted curve. Here we draw only 
the case of c = d = 1/2. One can see from the figure that, 
at 40 MeV, the neutron contribution to the total reaction 
probability is about two times of the proton contribution, 
while at 800 MeV the proton contribution exceeds that 
of the neutron. This reflects the behavior that the pn 
total cross section, 0^°*, is significantly larger than the 
pp total cross section, (7*°*, at low energy region. 

Figure [5] displays the reaction probability of proton- 
22 C reaction, similarly to that of proton- 12 C reaction. In 
contrast to the case of 12 C, the neutron contribution to 
the reaction probability is larger than that of the proton 
over all the energy range. Also, the reaction probabil- 
ity on the surface region comes mostly from the neutron 
contribution at all the energy range. This is due to the 



large extension of the neutron density, as shown in Fig. 
1. 

The difference between the reaction probability of 
proton- 12 C and that of proton- 22 C is as follows: Let the 
probing position of the proton be an impact parameter 
at which 2irbP(b) becomes a maximum. The probing po- 
sition for the case of proton- 12 C at 800 MeV is at 2.0 
fm, and the maximum height is about 10 fm, while they 
are about 2.7 fm and 13 fm in the case of 22 C. The reac- 
tion probability of proton- 12 C at 800 MeV reaches zero 
at about 6 fm, while in the case of proton- 22 C it reaches 
zero at about 10 fm. The major contribution comes from 
the region around the probing point, i.e., the surface. 

In order to compare the sensitivity of the proton and 
carbon probes to the nuclear surface, we plot, in Fig. [5[ 
ctr(&)/<tr for three nuclei, 12 C, 19 C and 22 C, of different 
features as a function of the impact parameter, b, where 
or (6) is defined similarly to Eq. (fT]) but the upper limit of 
the integration is limited to b. The nucleus 12 C is a stable 
nucleus which has almost the same proton and neutron 
distributions, 19 C is a good example of one-neutron halo 
nucleus, whereas 22 C is a two-neutron halo nucleus with 
a long neutron tail. 

As one can see from the figures, for each case, the 
major contribution comes from the surface region, which 
supports the above discussion. As a rough estimate of the 
extent to which the surface region is probed, we may use 
an impact parameter at which or (6) reaches 90% of or. 
Then we take the difference of such impact parameters, 
Ab, between 40 and 800 MeV incident energies. The 
increase of A6 for the change of incident energy from 
800 to 40 MeV is understood from the fact that the pn 
interaction becomes longer-ranged and stronger, which is 
reflected in the energy-dependence of (3 pn and 0^°* . 

First we focus on the reaction cross sections for the 
proton target. The A6 value increases from 0.6, 1.5 to 
1.9 fm as the neutron density becomes more widely dis- 
tributed for 12 C, 19 C and 22 C, respectively. This suggests 
that the proton target can probe the density distribution 
near the surface up to further distances as the interac- 
tion range increases. The corresponding Ab value for the 
12 C target case is 0.7, 1.3 and 1.6 fm for 12 C, 19 C and 
22 C, respectively. Comparing Ab values for the proton 
targets with those for 12 C targets, we can conclude that 
the 12 C target can probe the surface region equally to the 
proton target but is disadvantageous to probe the remote 
surface region of the spatially extended neutron distribu- 
tion, such as 22 C, compared to the proton. This is due 
to the fact that the proton and neutron distributions in 
12 C are very similar and that the nn(pp) interaction is 
shorter-ranged and weaker than the pn interaction. 

It would be possible to probe the outer region of the 
density distribution by the proton target especially at 
lower energy, but, at very low energy, we have to note 
that the long wavelength of the proton leads to a low 
resolution to the resultant density distributions, which 
may prevent us from studying minute structures of the 
outer density in detail. 
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FIG. 7: Reaction probability, Eq. (|16[) , times 2ivb for proton- 12 C. The solid curve is the total probability. The dashed curve is 
the neutron contribution, while the dotted curve is the proton contribution. The choice of c = d = 1/2 is made. 



VI. SUMMARY 

We have made a systematic analysis of the total reac- 
tion cross sections for the carbon isotopes incident on a 
proton target for wide energy range, in comparison with 
the features of a carbon target. We have predicted the 
reaction cross sections especially at 40 MeV where the 
experimental data have been measured at RIKEN. 

We have formulated this problem using the Glauber 
theory. The inputs are the parameters of nucleon-nucleon 
profile functions and the wave functions (densities) of 
carbon isotopes. The parameters of the nucleon-nucleon 
scattering are determined from the available experimen- 
tal data. The densities are generated using the wave 
functions of the Slater determinant, which we used in 
our previous work [14j . To go beyond that, we use a 
core+n model for odd nuclei and a core+n+n model for 
16 < 22 C nuclei. 

Having treated the interactions of proton-proton and 
proton-neutron separately, we have shown that the opti- 
cal limit approximation of the Glauber theory gives al- 



most the same results as the few-body calculation for the 
proton-nucleus reaction cross section over all the energy 
range used here. 

For 22 C, we generate several densities which are con- 
structed from the wave functions giving different sepa- 
ration energies of 0.489, 0.361, 0.232 and 0.122 MeV for 
the last two neutrons. All of them lie within the error 
bar of the experimental value, 0.423±1.140 MeV. The re- 
action cross sections calculated using these densities are 
957, 969, 985 and 1005 mb, respectively, for proton- 22 C 
at 40 MeV. Since the preliminary data of proton- 22 C re- 
action cross section has been reported to be around 1000 
mb with a large uncertainty, all of our predictions are 
consistent with the data, but the larger two values, 985 
and 1005 mb, would be favorable. If so, the data may 
suggest very small S^n- 

At around 100 MeV, the values of res defined by 
Eq. (|12p happen to agree reasonably well with r m listed 
in Table [Til which suggests that we may empirically ac- 
cess to the rms nuclear matter radii of carbon isotopes 
just by measuring or at 100 MeV. 



10 





FIG. 8: The same as in Fig. [7] but for proton- 22 C case. The choice of c = d = 1/2 is made. 



We have found a parameter free new relation, Eq. (fT5)) . 
It helps us to predict reaction cross sections for various 
isotopes at a given energy if the reaction cross section 
value of some stable isotope is available. 

Finally, we have made simple estimates for the contri- 
bution of the neutron and the proton to the total reac- 
tion cross sections. The major contribution to ctr comes 
from the surface region. Moreover, we have pointed out 
that a proton target can probe the surface region of the 
neutron-rich nuclei better than a 12 C target especially at 
lower incident energy. 
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Appendix: NN Scattering Amplitudes 

Here we discuss the parameterizations of nucleon- 
nucleon scattering amplitudes. In the text, we param- 
eterize it with a single Gaussian. We show here that 
the parameterization with double Gaussians gives numer- 
ically almost the same results for the reaction cross sec- 
tions as the single Gaussian, and validates our use of the 
single Gaussian prescription. 

We only show the case of pn scattering, because its con- 
tribution is more important for the neutron-rich isotopes 
than pp scattering, and also because, as we have discussed 
in Sec. V, the pn reaction dominates the proton-nucleus 
reaction cross sections especially at energies less than 100 
AMeV. 

In Fig. [TOl the numerical results of the pn elastic scat- 
tering differential cross sections calculated using the pa- 
rameters of Ref. [2(| are compared with the data. The 
numerical results of the single Gaussian are shown by 
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FIG. 9: Comparison of the calculated reaction cross section at a given impact parameter, b, in ratio to the total reaction cross 
section between proton- 12,19 ' 22 C and 12 c_ 12 . 19 > 22 q The solid curve is the calculation at 800 MeV, while the dashed curve is at 
40 MeV. The dotted line indicates 90% of or . 



the dashed curve. The agreement of the results with the 
data is reasonable, but not perfect especially in the for- 
ward direction. 

Although the agreement with the data of the pn elastic 
differential cross section is imperfect, the expression of 
the single Gaussian well reproduces the reaction cross 
sections as in Tabic fVTl 

For comparison, we perform fittings for the pn elastic 
scattering data in the energy range from 40 MeV to 100 
MeV, using the parameterization of double Gaussians. 
For this case, the profile function, r^y, for pp and pn 
scatterings, is parameterized in the form 

(18) 

where, a\, ai, Pi, a-i and /?2 are fitting parameters 
determined from the requirement: 1) The optical theo- 
rem is satisfied. 2) The ratio of the real to the imaginary 
part of the pn(pp) scattering amplitude in the forward 
direction reproduces the experimental values. 3) The to- 



tal elastic scattering cross section is equal to the total 
cross section. 4) The elastic scattering differential cross 
sections are reproduced. 

The fitting results using the double Gaussians are dis- 
played by the solid curves in Fig. |TD] The two sets of 
experimental data shown as 60 MeV in Fig. [10] are at 
62 MeV (open circle) and 63 MeV (thick dot). As for 
the differential cross sections, it seems that the double 
Gaussians give better results. 

In Table IVIl we compare the numerical results of 
proton- 12 C total reaction cross sections. The calcula- 
tions using the double Gaussians are the fit, while those 
using the single Gaussian are obtained by the use of the 
parameters of Ref. [2(| ■ The parameters of pp scattering 
are kept fixed. The experimental data are also shown in 
the table. The difference between the results of the reac- 
tion cross sections using these two parameterizations is a 
few % except at 60 MeV, but around this energy the data 
scatter widely and the difference between the numerical 
results and the data is not serious. 
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FIG. 10: The differential cross sections of the pn elastic scattering. The solid curves denote the fitting using double Gaussians, 
as explained in the text, while the dashed curves the results with the parameters listed in Table [I] The experimental data are 
taken from Ref. [281. 



Thus, we conclude that the parameters listed in Ta- 
ble Q] works fairly well for the reaction cross sections. For 
simplicity, we adopt the parameterization with the single 
Gaussian throughout this paper. 



E (MeV) Table [J present fit Exp. 



40 


432 


416 


371(11), 405(38) 


60 


359 


387 


310 (13) 


80 


314 


320 


279 (10) 


100 


284 


294 


275 (21) 



TABLE VI: Total reaction cross sections of proton- 1 C in mb 
calculated using the parameters given in Table [T] and that 
determined from two Gaussian fitting for pn elastic scattering 
data. The experimental data are taken from Refs. [26l. |27|. 
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